Astrophysical limits on light NMSSM neutralinos 
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It was recently shown that light LSP neutralinos could be found in the framework of the NMSSM. These 
candidates would escape known Particle Physics constraints even though they are relatively light. We now in- 
vestigate the astrophysical limits which can be set on these particles. We show, in particular, that the Fermi 
observation of dwarf spheroidal galaxies enable to constrain the parameter space associated with these candi- 
dates. Combined with the XENONIOO experimental limits, our results illustrate the complementarity between 
direct and indirect searches for dark matter. Yet, our findings also suggest that probing light neutralinos in the 
NMSSM scenario will be very difficult because the sensitivity of both dark matter direct and indirect detection 
experiments would have to be improved by at least six order of magnitude compared to present values in order 
to explore the entire parameter space. Finally, we show that the parameter space compatible with the CoGeNT 
signal (albeit disfavored by the XENONIOO limit) is not excluded by gamma nor radio observations. 



I. INTRODUCTION 

After several decades devoted to dark matter searches in 
underground, collider and spatial experiments, we are finally 
reaching an important cross road. Now that LHC is running, 
the Fermi and Planck experiments are delivering data and di- 
rect detection experiments have reached the level of sensitivity 
required to probe dark matter particles in the GeV-TeV range, 
we expect important developments in the dark matter field. 
In particular, it is likely that we obtain enough information 
in the next few years to determine whether supersymmetry 
manifests itself at the weak scale or not and, hence, whether 
neutralinos can constitute the dark matter. 

Until then, elucidating the nature of dark matter remains 
challenging. While ongoing experimental efforts try to close 
down the possible dark matter mass range, several dark matter 
direct detection experiments have announced events or signals 
which could point towards the existence of relatively light par- 
ticles CHSl. Although these seem also compatible with back- 
ground expectations, these claims have revived the theoretical 
interest for candidates in the GeV-10 GeV mass range and 
encouraged experiments to investigate the low energy range 
despite the lack of sensitivity at small recoil energies. 

From a theoretical point of view, motivating candidates in 
this mass range is not an easy task. For example, the low- 
est neutralino mass in the Minimal Supersymmetric Standard 
Model (MSSM) (with common slepton and squark masses and 
unification of the gluino and wino masses) now appears to 
be above 15 GeV |4, 5|. Relaxing the universality condition, 
scenarios where neutralinos could be below this value were 
found I16H9I, nevertheless the low mass range appears to be 
both statistically unlikely and challenged by direct detection 
and Higgs searches at colliders ^ . In view of these results, 
other candidates were investigated, in particular those in ex- 
tensions of the MSSM such as the neutralino in the Next- to- 
Minimal Supersymmetric Standard Model (NMSSM) ifTTI - 



^ Neutralinos lighter than 1 GeV have been shown to satisfy collider, astro- 
physical and cosmological constraints (To) but we do not consider these 
cases. 



[TtI, in the MSSM with an extended Higgs sector (TSl or the 
right-handed sneutrino in supersymmetric extensions |[T9ll2Ql . 

In a previous study, we demonstrated that NMSSM neu- 
tralinos could be very light |4|. In particular, we found 
that many points with a LSP mass between 1 and 15 GeV 
were actually compatible with known particle physics con- 
straints and possibly also with direct detection exclusion 
bounds. Here, we investigate whether these scenarios respect 
the most recent astrophysical limits since light (1-15 GeV) 
particles could overproduce the radio emission in the Milky 
Way (MW) and in galaxy clusters lEH |22J, gamma rays in 
dwarf spheroidal (dSph) galaxies [23l and antiprotons in the 
Milky Way (ni El. 

In Section II, we recall the parameter space associated with 
light NMSSM neutralinos (before applying the astrophysical 
and direct detection limits). We show, in particular, that some 
points have a very large annihilation cross section at small 
dark matter velocity despite a non negligible relic density. In 
Section III, we show how astrophysical and experimental lim- 
its cut into the parameter space, thus demonstrating the com- 
plementarity between direct and indirect detection searches. 
In Section IV, we estimate the radio flux expected for bench- 
mark points, selected so as to evade the Fermi and direct de- 
tection limits. We conclude in Section V. 



II. NEUTRALINO ANNIHILATION IN THE NMSSM 

In this section, we delineate the parameter space associ- 
ated with light NMSSM neutralinos and compute the dom- 
inant branching ratios. For this purpose, we consider ther- 
mal candidates and require that their energy density today 
is either equal to (or smaller than) the observed dark matter 
abundance, that is ^wmap^^ > ^%h^ > 10%^wmap^^ with 
^WMAp/^^ = 0.1131 ±0.0034 |25|. This basically constrains 
the total annihilation cross section of neutralinos in the pri- 
mordial Universe and rules out part of the parameter space. 
We also impose constraints from new particle searches at col- 
liders, from B -physics observables and from the muon anoma- 
lous magnetic moment as detailed in 1 4 1 . 

Once we obtained the configurations which satisfy all these 
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FIG. 1 : Relic density of light NMSSM neutralinos, the darker the dot 
the larger the likelihood. 



constraints, we can predict the energy spectrum of photons in 
the galaxy and dSph as well as the flux of cosmic rays in the 
MW. The comparison of the spin-independent cross section 
with the limits from CDMS |3 1 and XENONIOO (26] will en- 
able us to set additional constraints. In view of the tension at 
low dark matter mass between the results from the XENON, 
CDMS, CoGeNT |1| and DAMA/LIBRA |2| experiments, we 
will also investigate whether some scenarios can fall into the 
low mass and large cross section region or not. 

To perform this analysis, we have used the same 
Markov Chain Monte-Carlo code presented in |4| based 
on the NMSSMTools package (StI embedded in mi- 
crOMEGAs |28|. Dark matter observables were computed 
with micrOMEGAs 1291 [301. 



A. The model 

The Next-to-Minimal Supersymmetric Standard Model 
(NMSSM) is a simple extension of the MSSM that contains 
an additional gauge singlet superfield. The VeV of this sin- 
glet induces an effective ju term that is naturally of the or- 
der of the electroweak scale, thus providing a solution to 
the naturalness problem |31|. The model contains one ad- 
ditional neutralino state, the singlino, as well as three scalar 
{H\^H2-,H^) and two pseudoscalar (Ai,A2) Higgs bosons. An 
important feature of the model is that the singlet fields can 
be very light and escape the LEP bounds. This is because 
these fields mostly decouple from the SM fields |31|. This 
opens up the possibility for new annihilation mechanisms for 
light neutralinos, in particular through the exchange of light 
Higgs singlets as well as into light Higgs singlets f28|. The 
model that we consider has input parameters which are de- 
fined at the weak scale. The free parameters are taken to be 
the gaugino masses Mi ,M2 = M3/3, the Higgs sector param- 
eters ^, tanP, ^, k,A;^,Ak, a common mass for the sleptons mj 
and the squarks ruq as well as only one non-zero trilinear cou- 
pling. At, for more details see |i4J. We only consider scenarios 



with a neutralino LSP lighter than 15 GeV. 



B. Branching ratios in the Early Universe 

We start by computing the relic density for each candidate. 
The main assumption in these calculations is the conventional 
freeze-out mechanism. As shown, in Fig. [T] we have found 
many points with a relic density in the WMAP range. Hence, 
although the candidates that we are interested in are light, 
many have an acceptable relic density. In what follows, we 
do also include the points with a smaller relic density even 
though they can only partially contribute to the dark matter. 

We can now determine the Branching Ratios (BR) associ- 
ated with the different final states. The results are displayed in 
Fig.[2j The dominant annihilation channel is either into Higgs 
pairs, H\H\ and AiAi, or into Fermion pairs. We did not find 
any configurations with H\A\ in the final state as it would re- 
quire that both the scalar and pseudoscalar be very light. For 
Fermionic modes, the dominant channel is determined by the 
heaviest kinematically accessible Fermion. When the mass of 
the neutralino is smaller than 1 .7 GeV, the only possible final 
states are into light quarks ss and cc. If the neutralino mass 
is larger than 1.7 GeV but smaller than 4.2 GeV (m^), the 
dominant channel is XX at 90-100 %. Above the h mass, the 
dominant Fermionic final state is usually bb. However, the as- 
sociated branching ratio spans from below 1% to 100% as the 
Higgs mode can also contribute significantly. Hence, for neu- 
tralino masses above the b-quark mass, one expects also an- 
nihilations into H\H\^A\A\ as well as some contribution from 
XX and cc. 



C. Branching ratios in Milky Way and Dwarf Spheroidal 
galaxies 

We can now compute the annihilation cross section and 
branching ratios in the MW and dwarf spheroidal galaxies. 
In effect, this is equivalent to studying the impact of the dark 
matter velocity. To take into account the fact that neutralinos 
may not be responsible for all of the dark matter in the uni- 
verse and therefore that the halo may not be totally composed 
of neutralinos, we introduce the parameter 



^ = 1 if Q.j}? > ^WMAP/^^ 
= ^5^/z^/^wMAP^^ otherwise 



(1) 



where ^/^wmap corresponds to the lower value of the WMAP 
measurement. In what follows we will always rescale the local 
neutralino density by the factor, p^^^ = ^Pdm- 

The total annihilation cross section spans several orders of 
magnitude as displayed in Fig. [3] In some cases, it can be 
strongly enhanced with respect to its value in the Primordial 
Universe. This "boost" can occur when the annihilation pro- 
ceeds through a s-channel exchange of a pseudoscalar Higgs 
particle near resonance, the cross section is then sensitive to 
the thermal kinetic energy : at small velocities, one gets the 
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FIG. 2: Branching ratios associated with neutrahno pair annihila- 
tions in the Early Universe. The color code is associated with the 
constraints from dSph (red) and direct detection (yellow), see sec- 
tion III. 



full resonance enhancement while at v ~ c, one only catches 
the tail of the resonance 1321 [33ll . 

The associated cross section is proportionnal to 



va(v) 



1 



1 



16m4 (v2/4 + A)2 + r2(l-A)/4m2 



(2) 



where A = 1 — m^/4m^. At v ^ 0, it is strongly enhanced as 
compared to its value at freeze-out for A, Fa ^ 1. 

To give a more quantitative estimate of this effect, let 
us consider one allowed scenario with m^^ = 10.08 GeV, 
niA^ = 20.12 GeV, T = 1.1 x 10""^ GeV and compute the ratio 
of the thermally averaged cross section at a given tempera- 
ture, (av)(r) to the value at a typical freeze-out temperature, 
(av) (r = m/20). The enhancement factor reaches two orders 
of magnitude and depends mostly on A since the term in Fa in 
Eq.[2]is negligible. A small variation in A can lead to an even 
larger enhancement factor, see Fig. |4] where the enhancement 




FIG. 3: Rescaled neutralino annihilation cross section in the galaxy 
versus the neutralino mass, the points which overpredict the gamma 
ray flux in dSph are in red, see section III. 



factor is displayed for different values of A. 
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FIG. 4: Ratio of (av) to {ov){T = m/20) as a function of the tem- 
perature. 

As a result, we find that the neutralino annihilation branch- 
ing ratios in the MW or dSph differ from those found in the 
early Universe. Our results are summarized in Fig. [5] where 
we only display the final states which opened up at low veloc- 
ities. The AiAi and //i//i are no longer possible final states 
because the associated cross sections are both suppressed. On 
the other hand, the branching ratios for the cc and XX now 
reach unity for many points where m^ > rric and m^ > m^ re- 
spectively. 

Many scenarios have indeed a very small annihilation cross 
section in dwarf galaxies due to a p-wave suppression factor. 
This affects, in particular, the annihilation processes which 
were dominated in the Early Universe by a Hi resonance de- 
caying into Fermion pairs as well as annihilation into light 
Higgs final states through the t and u channel neutralino ex- 
change. For these configurations, channels such as the t- 



channel sFermion exchange, which were subdominant in the 
Early Universe, become important when v ^ 0. 



hilation into all SM final states, is thus given by 




FIG. 5: Branching ratios associated with neutralino pair annihilations 
at low velocities. 



In principle, we can make generic predictions for all the 
configurations with a BR = 1 for a given channel since the 
annihilation cross section can be factorized. However, most 
of the dark matter scenarios that we have found have a mix- 
ture of final states. Hence, in the following, we will compute 
individually the gamma and cosmic ray fluxes for each point 
considered. 



III. ASTROPHYSICAL LIMITS 

We can now investigate the astrophysical limits which ap- 
ply to light neutralinos. 



A. Gamma rays in dSph 

Objects such as dSph are meant to be dark matter domi- 
nated and represent therefore in principle a good target for 
dark matter searches (231 [34, 35 1. Dark matter annihilations 
are expected to produce quarks and/or taus in the final state 
which, in turn after hadronization, produce gamma rays. Pho- 
tons can also be produced directly as radiation from an in- 
ternal line or from a final state. Finally it is also possible 
that neutralino annihilate directly into photon pairs through 
a loop-induced process, but the cross sections are typically 
small |36|. 

The flux of gamma rays originating from dark matter anni- 
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dNi 
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^^ jdm^l^diyf)) 



(3) 



where {(5v)tot is the total annihilation cross section and BRi is 
the fraction into a given SM final state /, pDM is the dark matter 
energy density, dNi/dE is the number of photons produced 
after hadronization (or radiation even though this process is 
generally subdominant) in terms of the energy E, /(\|/) is the 
line of sight in the \|/ direction. The factor 1/2 accounts for 
the Majorana nature of the neutralino. 

This flux has to be integrated over the Fermi angular reso- 
lution. For an angular region of diameter 0.5°, as considered 
by the Fermi experiment |37|, one is not sensitive to the inner 
slope of the dark matter halo. Hence we can safely consider a 
NFW dark matter halo profile. The value of the integral along 
the line of sight and averaged over the resolution is taken from 
Table 4 in [37 1. The energy dependent part of the differential 
flux is then computed with micrOMEGAs and integrated from 
0.1GeV<£<m,o. 

M 

For each point found by our MCMC, we have computed 
the gamma ray flux expected in the eight dwarfs considered 
by the Fermi experiment. We then compare this value with 
the Fermi-LAT 95% limits |37l. Our results are displayed 
in Fig. [6] for Draco which gives the most sensitive limits. As 
expected, excluded models are those with the largest annihi- 
lation cross section (irrespective of the dominant annihilation 
channel, see also Fig. |3]) and for which the mass difference 
between the LSP and the pseudoscalar involves a fair amount 
of fine tuning. 

Note that the exclusion limit varies only slightly with 
the mass. Basically, it corresponds to i^ov/m^ > 2 — 5 x 
10~^^cm^s~^GeV~^. Also the criterion for exclusion is based 
on a comparison of the computed flux with the Draco limit 
taking into account the la error bars in the integral over the 
DM density distribution. Thus, it is a conservative exclusion 
criterion. However, many points excluded by DRACO are 
also excluded by other astrophysical observations (including 
by the flux of antiprotons measured at Earth position). The 
implications for the Higgs spectrum will be discussed in sub- 
section |IIlDl 



B. Antiproton and positron fluxes in the Milky Way 

We can now investigate the cosmic ray production in the 
MW. Calculations are similar to gamma rays except that now 
we focus on positron and antiproton production. We also ex- 
pect that the best signals are associated with the largest anni- 
hilation cross section. 

Since the limits from dSphs are already very constraining, 
we shall consider other astrophysical bounds as a complemen- 
tary tool and thus apply them only to the scenarios which sur- 
vived the Fermi constraints. For this purpose, we group all 
the scenarios in bins of 1 GeV for neutralino masses ranging 
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FIG. 6: Predicted photon flux as a function of the LSP mass. The 
horizontal Unes correspond to the Fermi Umit including the la error 
bars in the integral over the DM density distribution. Points excluded 
by XENONIOO are in yellow. 



FIG. 7: Cosmic ray predictions for the antiproton spectra for 13 
benchmark points in Table [l| The first point in this table is ignored 
as it is too light to give antiprotons. The maximal energy for each 
spectrum corresponds to the neutralino mass. The background from 
secondaries |38 1 (dash) is also displayed. 



from 1 to 15 GeV and, for each of the 14 bins, we select the 
point with the maximum value for ^^av/m^o ^^^ with a good 

likelihood (defined as 2 > 032Qmax) and which, yet, is safe 
with respect to the Draco limits. This leads to 14 benchmark 
points. The dominant annihilation channel for all of them is 
into bb for m^o > 4.2GeV and into XX otherwise. The bench- 

M 

marks are listed in Table HI 

The results for the antiproton fluxes are displayed in Fig.|7] 
and compared to the background spectrum taken from the ap- 
proximate analytical formulae in Ref 1,38.1 . This background 
by itself provides a good fit to the data. The fluxes computed 
are close to one order of magnitude below the background 
from secondaries at energies below ^ 2 GeV and drop rapidly 
at higher energies. We therefore conclude that one cannot fur- 
ther constrain these scenarios by measuring the antiproton flux 
as the uncertainty in the background calculation (which ex- 
ceeds 10%) is always larger than the signal. 

Note that to compute these fluxes we have set the propaga- 
tion parameters to the default values in micrOMEGAs, i.e. the 
MED set of parameters (see Ref. 1391 ). For a different set of 
propagation parameters, the expected fluxes can increase. In- 
deed we find that the fluxes reach at most the background level 
with the MAX set of propagation parameters (corresponding 
to a larger diffusive zone, see Ref. 1 39", "401). 

For completeness, we have also computed the positron 
fluxes for the 14 benchmark points, using the MED propaga- 
tion parameter set |40|. The fluxes are always at least two or- 
ders of magnitude below the background |41 1 , see Fig. [8] This 
was to be expected since the scenarios that we have selected 
predict dominant quark final states; hence a better signature in 
antiprotons than in positrons. Note that the XX channel, which 
leads to a hard positron spectrum, is dominant either at very 
low masses or when the cross section is very small. 
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FIG. 8: Predictions for the positron spectra as compared to the back- 
ground |41 1 (dash) for 14 benchmark points in Table [l| The maximal 
energy for each spectrum corresponds to the neutralino mass. 



C. Comparison with direct detection 

In the light neutralino scenarios there is a good comple- 
mentarity between gamma ray searches and direct searches. 
Indeed many scenarios which predict a spin independent cross 
section below the XENONIOO exclusion curve overpredict the 
gamma ray flux in dSph galaxies, see Fig. [9] This is also il- 
lustrated in the correlation plot displayed in upper panel of 
Fig. 10 where we show the gamma ray flux as a function of 



the spin independent cross section. Clearly, the Fermi dSph 
limits constrain scenarios where the spin-independent cross 



niyO 

[GeV] 


^ 


(av) X 102'7 
[cm^s-i] 


BR,, 


BRbb 


BRss 


R 


0.976 


0.373 


0.209 








0.997 





2.409 


1.00 


0.297 


0.964 





0.026 


0.040 


3.342 


0.935 


0.345 


0.972 





0.018 


0.044 


4.885 


0.465 


3.298 


0.0970 


0.901 


0.0016 


0.041 


5.626 


0.376 


5.389 


0.0698 


0.929 


0.0011 


0.040 


6.551 


0.528 


3.547 


0.0618 


0.937 





0.046 


7.101 


0.689 


2.425 


0.0586 


0.940 





0.050 


8.513 


0.829 


2.161 


0.0416 


0.958 





0.055 


9.274 


0.827 


2.497 


0.0533 


0.946 





0.060 


10.27 


0.906 


2.323 


0.0634 


0.935 





0.063 


11.50 


0.960 


2.575 


0.0611 


0.937 





0.074 


12.74 


0.955 


3.224 


0.102 


0.897 





0.088 


13.51 


0.558 


9.571 


0.0781 


0.921 





0.085 


14.48 


0.147 


148.4 


0.0748 


0.924 





0.088 



TABLE I: Benchmark points: main characteristics and ratio of the 
dark radio emissivity at 330MHz to observation (R). 
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FIG. 10: Correlation between the gamma ray flux and spin inde- 
pendent cross section of NMSSM neutralinos. Top: all points are 
included. The yellow points correspond to scenarios with a too large 
spin independent cross section. The dashed line corresponds to the 
Fermi limit for the Draco dSph. Bottom: all points overpredicting 
the gamma ray flux or with a too large spin independent cross sec- 
tion and which do not completely explain the dark matter today have 
been removed. 



FIG. 9: Spin independent cross section versus the neutralino mass. 
In red are the points which over predict the gamma ray flux in dSph. 



section is smaller than the latest XENON 100 limit while the 
XENON 100 limits exclude points where the gamma ray flux 
in dSph is not yet accessible by the Fermi searches. In the 
framework of the NMSSM, this complementarity is directly 
connected to the light Higgs spectrum as discussed in the next 
subsection. 

If we now remove the points which do not have the cor- 
rect abundance today and exclude the points which pro- 
duce too many gamma rays in Draco and non observed 



events in XENON 100 (see Fig[T0| lower panel), we ob- 
tain that, statistically, light neutralinos are likely to produce 
[10~^^, 10~^^] y/cm^/s and have a spin independent cross sec- 
tion of [IQ-^^, IQ-^^] cm^. 



We have also computed the gamma ray flux for points 
which are in the region favoured by the CoGeNT experiment 
IT]. All these points lie in the region excluded by XenonlOO. 
However, since CoGeNT claims detection at 2 a of an annual 
modulation signal ll42ll . it is worth investigating the astrophys- 
ical limits for such candidates. 

Since we have demonstrated that indirect and direct detec- 
tion experiments were probing different regions of the param- 
eter space and these candidates are within XENON 100 sen- 
sitivity, we do not expect that they produce large gamma ray 
and cosmic ray fluxes. However, to check this statement, we 
shall consider three benchmark points (cf Table lll|. 



For these points, we found ^ Gv/n 



< 6 X 



10~ cm s~ GeV~ which is one or two orders of magnitude 
below the Draco limit in section |III A| Hence, it seems that 
NMSSM neutralinos in the CoGeNT region are not excluded 
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TABLE II: Three examples of NMSSM points falling in the CoGenT 
contour 1 1 1 in the ^a^^ vs. rUyQ plane. For all of them we have set 
M3 = 3 X M2 and Aij=A^ = 0. All quantities are expressed in GeV 
units. 



by the indirect detection limits for the moment. 

This is an important point nevertheless. Indeed, it will be 
challenging for both dark matter direct and indirect detection 
experiments to reach the level of sensitivity which is required 
to completely probe the region of the parameter space with 
light neutralinos. In addition, new sources of background 
(such as the neutrino background for the direct detection ex- 
periments) may also weaken the analysis. 



D. Implication for particle physics 

In the previous subsections, we have demonstrated that 
the Fermi dSph limits were setting stringent limits on the 
NMSSM parameter space and were complementary to dark 
matter direct detection searches. We can now examine the 
impact of these limits on the Higgs sector and on B -physics 
observables. 
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FIG. 11: Correlation in the (mAy.mHi) plane. All the safe configura- 
tions (green) as well as the excluded points by XENON 100 (yellow) 
and dSph (red) are displayed. 



Efficient neutralino annihilation in the Early Universe re- 
quires at least one light Higgs (m//^ ,mAi < 30 GeV for m^ < 
15 GeV) as illustrated in Fig. [TT] Astrophysical limits then 
apply in two distinct regions of the ma^ — mn^ plane. The 
first region corresponds to a light Hi (niHi ^ i^^^^] GeV) and 



to heavier A 1 (with jua^ G [10, 1000] GeV). In this region the 
spin-independent cross section can become very large, which 
is in conflict with XENONIOO data (in yellow in Fig.[TT]). In- 
deed, as can be seen in Fig. 12 larger spin independent cross 
section are found for light Hi . This is because the scalar ex- 
change contribution to the cross section goes as 1/m^^ . Note 
that because sufficiently large couplings of the light Higgs to 
the LSP and to quarks in the nucleon are necessary to have 
a large SI cross section IH, many points with light Hi are 
not excluded. The second region corresponds to a relatively 
light Ai (rriA, G [10,30] GeV) and rriH, G [20, 100] GeV Here 
the neutralino pair annihilations (which proceed through the 
exchange of an A 1 in the s-channel) become singular when 
vjm -^ and 1 — mAi/2m^ << 1 and thus can produce too 
many gamma rays in dwarf galaxies (see points in red in 



Fig. 12). 



E 
o 



CO 

D 



^0-48 



■ ■• .■'■-■*i 




I I I 

10 



10^ 

MnJGeV] 



FIG. 12: Spin independent cross section in terms of m//^ , same color 
code as Fig.[TT] 



The observables in the B- sector have been used to constrain 
the parameter space. The LHCb experiment, which is now 
taking data, will measure these observables with increased 
precision, it is therefore interesting to examine whether this 
will probe further our scenarios. For example the branch- 
ing ratio for Bs -^ jn^jn' is expected to be powerful in prob- 
ing scenarios with a light doublet Higgs and large values of 
tan p. Only a fraction of our scenarios have such charac- 
teristics. The predictions for B{Bs -^ !^^l^~) are displayed 



in Fig. 13 together with the expected limit of LHCb with 
L = Ifh^ |43 j. A signal is expected only for a small fraction 
of the scenarios while in some cases the predictions are sup- 
pressed as compared to the SM expectation {B{Bs -^ A^^A/~) = 
3.6 ±0.4 X 10~^) EH. This implies that light neutralino sce- 
narios cannot be probed entirely with this observable. 
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FIG. 13: Predictions for the branching ratio B{Bs -^ IJ^'^IJ~) in terms 
of rrij^^ together with the Tevatron Umit and the LHCb projected hmit, 
same color code as Fig.pT] 



IV. RADIO EMISSION IN THE GALACTIC CENTRE AND 
THE COMA CLUSTER 

Since light particles eventually produce electrons in the 
MW, one expects a significant radio emission in the galactic 
centre as well as in galaxy clusters as first pointed out in 1211 
and discussed further in e.g. Ref. ll22ll45HJ8l . 

This is true, in particular, for the points with large cross sec- 
tions at V ^ and ^ ^ 1 (i.e. canonical values in the primor- 
dial Universe). Since these points have been already excluded 
by using the Fermi dSph limits, we shall rather concentrate 
on the part of the parameter space that is left after applying 
the gamma ray constraint. We thus compute the radio flux 
expected for the benchmark points displayed in Table [l| 

We repeat the same calculations as in Ref. |22| where we 
assume the MED set of propagation parameters and set the 
magnetic field, B, to 20 /iG. This value is slightly higher than 
that derived in [49 1 for very small scales but it is still con- 
servative enough with respect to the very large value gener- 
ally considered for the magnetic field in the galactic centre. 
Choosing the MAX set of propagation parameter would in 
fact decrease the intensity of the radio emission expected in 
the galactic centre but it would also lead to a broader radio 
emission in the galaxy. On the contrary, for the MIN set of 
propagation parameters, one expects a brighter emission in 
the galactic centre which should be easier to constrain in prin- 
ciple. 

We focus on one radio frequency, namely 330 MHz where 
the emission is about 360 Jy (for an angular resolution of 7') 
and compute the ratio R of the "dark" emissivity to the obser- 
vation at this frequency in the galactic centre, for each of the 
points mentioned above. The results are displayed in Table |lj 
in the last column. We find that all the benchmark points have 
a negligible radio flux. This can be explained by the fact that, 
in all these scenarios, the total number of electron produced 



in the 0.5 to 15 GeV range is about unity (a minimal injection 
energy of 1 GeV is required for B = 20 juG to produce syn- 
chrotron emission at 330 MHz). Indeed, all the benchmark 
points annihilate preferably into b or x-pairs. This leads to a 
very small number of electrons per energy unit for £^ > 1 GeV. 

Note that the determination of the magnetic field value is 
of particular importance because it changes the synchrotron 
emission (hence the radio predictions) and the losses. How- 
ever, it also determines other possible signatures such as an 
anomalous submillimetre emission in the galaxy which could 
contribute as a new source of foreground in Cosmological Mi- 
crowave Background studies. For example, to get a signal at 
33 GHz (relevant for both WMAP and Planck) from a 10 GeV 
dark matter particle, the magnetic field should be greater than 
25^G |50|. If such a value is indeed attained in the galactic 
centre or in the outer parts, one can hope to correlate "dark" 
radio emission with WMAP and Planck observations e.g. 1501 - 
"531. 

One can also estimate the radio emission in clusters of 
galaxies to determine whether light dark matter particles are 
potentially constrained by radio observations in these objects 
(21 , 22 , 45 1 . To compute the flux, we consider the Coma clus- 
ter and extend the procedure described in (21] |22l to account 
for the energy distribution of electrons. We assume a NFW 
profile |54| with po = 4.4 10"^ GeV/cm^ r, = 400 kpc, a de- 
tector angular resolution of 1 deg, a magnetic field of 4.7 juG 
and a density of electrons of 3 10~^ cm~^. 

The results for our benchmark points are displayed in 



Fig 



14 



^. Our prediction shows that none of these points are 
excluded by the observation in the Coma cluster, as expected 
already from our computations in the galaxy. 

These results do not account for substructures. Also, they 
assume a specific value of the magnetic field. If we increase 
this value by a significant amount, the radio flux becomes 
larger and the emission becomes also possible at higher fre- 
quencies. For example, increasing the value of the magnetic 
field up to 12 /i G increases the radio flux at 4.58 GHz by a 
factor ~ 4.5 for the last benchmark point (corresponding to a 
candidate with 14.48 GeV mass). Still this is not enough to 
rule out this candidate. 



V. CONCLUSION 

In this article, we have investigated the astrophysical lim- 
its on light NMSSM neutralinos. We have shown that the 
Fermi searches for gamma rays from dwarf Spheroidal galax- 
ies set stringent limits on the NMSSM parameter space. In 
addition, combining these results with indirect and direct de- 
tection searches restrict the parameter space even further. Yet 
light neutralinos are not ruled out. 

It may be extremely difficult to completely probe light neu- 
tralinos in the forthcoming future. Indeed, this would require 



^ We have checked that we could recover the results in |451 for a 40 GeV 
candidates annihilating into a.bb pair and with the same cross section. 
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and investigate distinctive signatures involving light neutrali- 
nos, for example invisible decay modes of the Higgs 1 13 1. Fur- 
thermore the mass of a light neutralino, produced in the decay 
of a selectron, could be determined at a future linear collider, 
for this selectrons have to be light enough to be pair produced 
1561 . Light DM particles can also have an impact on the CMB 
ll2n[57lJ59l. In Ref 158", '591, it was shown that their annihila- 
tion cross section is constrained by WMAP to be around the 
usual freeze-out cross section unless DM particles annihilate 
mainly into neutrinos, as in |60| for even lighter DM candi- 
dates. However, the Planck satellite should further probe these 
scenarios, including by detecting an anomalous synchrotron 
emission |50|. Nonetheless, at present, our results demon- 
strate that the complementarity between direct and indirect 
searches is a reality and a necessity, at least in the framework 
of supersymmetric scenarios. 



FIG. 14: Radio flux in the Coma cluster for 13 benchmark points 
displayed in Table. |l| from the second lightest (left) to the heaviest 
(right) candidates. The dataset are taken from the references in 1.55 J . 



to improve both direct and indirect detection experiments sen- 
sitivity by several order of magnitudes while experiments may 
run into new background sources. 

The LHC could nevertheless provide crucial information on 
the existence and on the spectrum of supersymmetric particles 
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